Measurement of the t] mass at KLOE. 
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An integrated luminosity of 410 pb~^, corresponding to ~ 17 million of rj events, has been analyzed to measure 
the rj mass using the decay 77 — » 77. The measurement is insensitive to the calorimeter energy calibration and the 
systematic error on the measurement is dominated by the uniformity of the detector response. As a cross check 
of the method the tt" mass from the decay 4> 7r°7, tt'^ 77 has been measured and it is in agreement with the 
most accurate previous determinations. 

The result obtained is = 547.873 ± O.OOTatat. ±0.031ayst, MeV, that is today best measurement of the rj mass. 



1. Introduction 

The KLOE experiment is performed at the 
Frascati </> factory DA$NE[T]. DA$NE is a high 
luminosity e"'",e~ colhder working at -y/s ~ 1020 
MeV, corresponding to the cf) meson mass. In the 
whole period of data taking (2001 - 2006) KLOE 
has collected an integrated luminosity of 2.5 fb~^, 
corresponding to about 8 billions of 4> produced 
and 100 millions of 77 mesons through the elec- 
tromagnetic decay (j) rij. 

The KLOE detector consists of a large cylin- 
drical drift chamber [5], DC, surrounded by a 
lead/scintillating- fiber sampling calorimeter [3], 
EMC, both immersed in a solenoidal magnetic 
field of 0.52 T with the axis parallel to the beams, 
z in the following. The DC momentum resolu- 
tion for charged particles is Sp±/p±=OA%. The 
calorimeter is divided into a barrel and two end- 
caps, and covers 98% of the total solid angle. 
Photon energies and arrival times are measured 
with resolutions as/E = 0.057/ (GeV) and 
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at = 54 ps/y^E (GeV) © 50 ps, respectively. 
Photon-shower centroid positions are measured 
with an accuracy of a = 1 cm/ ^JE (GeV) along 
the fibers, and 1 cm in the transverse direction. 
A photon is defined as a cluster of energy deposits 
in the calorimeter elements that is not associated 
to a charged particle. We require the distance 
between the cluster centroid and the nearest en- 
try point of extrapolated tracks be greater than 
3x(t(z, 0), where 4> is the azimuthal angle. 

The trigger [3] uses information from both the 
calorimeter and the drift chamber. The EMC 
trigger requires two local energy deposits above 
threshold {E > 50 MeV in the barrel, E > 150 
MeV in the endcaps). The trigger has a large 
time spread with respect to the time distance be- 
tween consecutive beam crossings. It is however 
synchronized with the machine radio frequency 
divided by four, Tsync=10.85 ns, with an accu- 
racy of 50 ps. For the 2001-2002 data taking, the 
bunch crossing period was T=5.43 ns. The time 
(To) of the bunch crossing producing an event is 
determined offline during event reconstruction. 
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2. Measurement of the 77 mass. 

The value of the 77 mass has been recently mea- 
sured with high precision by two collaborations 
NA48[S] {rrir, = 547.843 ± 0.030 ± 0.041 MeV/c^) 
and GEM|] (to^ = 547.311 ± 0.028 ± 0.032 
MeV/c^) using different techniques and produc- 
tion reactions. The two measurements differ by 
more than eight standard deviations from each 
other. The GEM measurement is in agreement 
with the older ones [7] while the NA48 measure- 
ment is higher. For this reason it is interesting 
to provide a further measurement of comparable 
precision in order to clarify the experimental sit- 
uation. 

3. Measurement method. 

We measure the mass studying the decay 
777, 77 — + 77. A kinematic fit is performed im- 
posing the 4 constraints given by the energy- 
momentum conservation. Sing there are three 
photons and there are 4 constraints, the fit over- 
constrains the energies of the photons that are, 
practically, determined by the position of the 
clusters in the calorimeter. The inputs of the 
fit are the energy, the position and the time of 
the calorimeter clusters, the mean position of the 
e+e" interaction point, the total four-momentum 
of the colliding e+e~ pair. Each of these variables 
is determined run by run using e+e^ e+e^ 
events ( almost 90000 events for each run, allow- 
ing a very precise determination of the relevant 
parameters). 

4. Selection. 

The — > ?77 events are selected by requiring 
at least three energy deposits in the barrel with 
a polar angle 9j : 50° < 6^ < 130°, not associ- 
ated to a charged track. A kinematic fit imposing 
energy-momentum conservation and time of flight 
of photons equal to the velocity of light is done for 
all 3 7's combination of N detected photons. The 
combination with the lowest is chosen as a can- 
didate event if x'^ < 35. The events surviving the 
cuts are shown in figH] where the Dalitz plot is 
shown. Three bands are clearly visible. The band 
at low 777,^ is given by the </> — > tt'^j, tt" — > 77, 
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Figure 1. Dalitz plot distribution of the selected 
37 events. The photons are sorted according their 
energies Ei < E2 < E3. 



while the other two bands are (j) ~^ VI ~^ 77 
events. 

With the cut shown in the Dalitz plot we select 
a pure sample of rj, tt*' —> 77 events. The resulting 
spectrum (figEl) can be fitted well with a 
single gaussian of ct '--^ 2.1MeV/c^ . 

In order to determine the systematic error we 
have evaluated the uncertainities on all the quan- 
tities used in the kinematic fit and the effect on 
the fitted value. 

A sample of e'^e^ — )■ 7r+7r~7 events has been 
used to check the mean position of the interac- 
tion point, the energy response of the calorime- 
ter and the alignment of the calorimeter with the 
Drift Chamber. The mean position of the inter- 
action point, determined run by run using the 
e~^e~ — > e'^e~ events, has been compared with 
the reconstructed tt^tt" vertex. The difference 
between this two values has been computed run 
by run and the spread of the distribution is used 
for the systematic error on the determination of 
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Figure 2. Invariant mass distribution of the 7172 
pair. The fitted function is the continuous hue. 



the LP mean position. Since the LP position is 
determined with the tracks reconstructed in the 
Drift Chamber, small displacement between the 
DC and the calorimeter can affect position mea- 
surement of this point in the calorimeter refer- 
ence frame. To check disalignment between the 
calorimeter and the Drift Chamber, the tt^ and 
TT"*" tracks of the 7r"'",7r~7 events were extrapo- 
lated to the calorimeter and the closest approach 
point to the cluster centroid was determined. The 
difference in the position Xdu — Xcst were deter- 
mined and the spread of these values are taken 
as systematic error on DC-Calo allignment. A 
small correction of 1.1 mm in the Calorimeter 
position along the y direction, the vertical, and 
2 mm along the z direction (the direction of the 
beam axis) was applied. 

The absolute energy scale of the calorimeter 
and the linearity of the energy response was 
checked using the e^e~ e~^e~j events and the 
TT+Tf-'-y events. The energy of the 7 can be deter- 
mined using the two charged tracks in the Drift 



Chamber and then compared to the reconstructed 
cluster energy. A linearity better than 2 % was 
found while the absolute scale was found to be 
calibrated at better than 1 %. These systematic 
uncertainties result in just 4 keV for the scale 
and 4 keV for the linearity on the value of the re- 
constructed mass. The measurement shows very 
small sensitivity to the calorimeter calibration be- 
cause, as explained before, the kinematic fit over- 
constraints the photon energy with the cluster po- 
sitions. 

For this reason it is important to evaluate the sys- 
tematic error due to the misalignment of single 
modules in the calorimeter. This has been done 
evaluating the value of the mass as a functon of 
the position of the photons in the calorimeter. A 
spread of about 10-15 keV was found and assumed 
as systematic error. 

Systematics due to the particular choice of the 
cut on the Dalitz plot shown in fig[l]was also de- 
termined to be 12 keV, while the cut on the 
pratically doesn't have any effect on the value of 
the 77 mass (0.7 keV). 

The measured value of the mass is, instead, very 
sensitive to the energy in the center of mass of the 
r]"f system. Due to initial state radiation emission 
(ISR) the available center of mass energy is a bit 
lower than the ^/s of the e^e~ beams measured 
using e+e^ — s- e+e^ events. A variation of 100 
keV of the measured mass value is predicted by 
the MC simulation. Since this correction is re- 
altively large we have checked the simulation of 
ISR emission in the MC in the following way: 

• the correction to apply to the fitted value in 
order to obtain the real value of the mass 
has been determined as a function of i/s 
and shown in figO 

• the whole data taking has been devided in 
ranges of y/s using 8 points for the on peak 
data and two off-peak points at 1017 MeV 
and 1022 MeV; 

• the value of the mass obtained for each 
value of ^/s has been corrected according 
the MC prediction; 

• the residual spread of these points has been 
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Figure 3. correction to the rj mass evaluated by 
MC as a function of -^s. Points without error 
bars: MC; continuous hue: fitted curve to the 
MC points; points with error bars: DATA. 
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Figure 4. as a function of ^/s after the ISR 
correction. The middle line is the mean value 
while the two line above and below are +/— 1 a 
systematic. 



taken as systematic error (8 keV) (see fig. 

El). 



All these studies have been done also for the tt" 
mass using the cj) — s- tt'^j events. The ratio of the 
two masses r — r7i,,/m^o has also been studied. 
All the contributions to the systematic error are 
summarized in table [TJ 

Finally stability versus run conditions are 
checked dividing the whole period of DATA tak- 
ing 2001-2002 in 8 sub-periods and determining 
the value of 771,^0, to,, and rrirj/rn^o as a function 
of the acquisition period. The results are shown 
in fig. El 

The values shown in the figure have been fit- 
ted with a constant and the results of the fits are 
shown in tab. |21 



Fit results 




Value 


Error 


X^/n.d.f 


C.L 


nrin 


547791 kcV 


7 keV 


6.9/7 


45% 




134886 keV 


12 keV 


7.7/7 


34% 




4.0610 


0.0004 


8.9/7 


26% 



Table 2 

Fit results to the 8 data points. 
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systematic effect 


(keV) 


TO^o (keV) 


mri/m^a X 10 ^ 


Calorimeter energy scale 


4 


1 


5.6 


Calorimeter not linearity 


4 


11 


31 


Vertex position 


4 


6 


19 


Angular uniformity (j) 


15 


12 


37 


Angular uniformity 6 


10 


44 


120 


ISR effect 


8 


9 


28 


Dalitz plot cut slope 


12 


4 


15 


Dalitz plot cut constant 


12 


1.9 


10 


cut 


0.7 


4 


13 


overall 


27 


49 


136 



Table 1 

Summary of all evaluated systematic effects. The Dalitz plot slope and constant refers to the change of 
slope or intercept of the stright line cut in fig. [1] 



5. Computation of the final result. 

We have two different ways to extract the value 
of the rj mass. We can use the ratio m^/m^o 
obtaining: 

^ = 4.0610 ± G.0004(stat.) ± 0.0014(syst.) 

from which using the PDG2006 f7] value of the 
7r° mass [mT^o — 134976.6 ± 0.6 keV) we obtain: 
m„ = 548140 ± 50(,tat.) ± 190(,y,t.) keV. 
Alternatively we can use directly the value of to^ 
coming from the fit. For this purpose we need 
to calibrate the ^/s with high precision (771^/771,^0 
is in fact -/s independent while the two values of 

and m^r" a-re fully correlated with the -^s mea- 
suremnt). The absolute ^/s scale has been cali- 
brated using the value measured by CMD-2 
[8j. To this pourpose the cross section e+e^ 
(j) — > KsKl has been measured as a function of 
the -v/s using the two off-peak points at ^/s — 
IQlTMeV and ^/s = 1022 MeV together with the 
on-peak data. The (j> resonance curve has been fit- 
ted using the CMD-2 parametrization that takes 
in account ISR and threshold effect in KsKl pro- 
duction. The central value of the 4> mass has been 
measured obtaining = 1019.329 ± 0.011 MeV 
[9]. The difference respect to the CMD-2 value 

CMD-2 = 1019.483 ± O.Oll.ta*. ± 0.025,^,,*. 
sets our absolute -^/i calibration. 

This means that our measurements can 
be regarded as a measurement of the ratio 



Tnn/m^CMD2 and m^o/m^cMD2'- 

m^/m4,cMD2 = 0.537403 ± 0.000007 (stat.) 

±0.000026 (syst.) ± 0.000006 (m^, stat.) 

and 

m^o/m^cMD2 = 0.132328 ± 0.000012 (stat.) 

±0.000048 (syst.) ± 0.000001 (m^stat.) 

We can extract the 7r° and 77 mass using the value 
of CMD-2 777,0 obtaining: 

= 134.906 ±0.012(stat.)±0.048(syst.) MeV 
771^ = 547.873 ±0.007(stat.)±0.031(syst.) MeV 

where the statistical and sytematic error have 
been independently summed in quadrature. 

The tt" mass is in agreement at 1.4 cr level with 
the PDG06 value, validating the whole procedure. 

The value of the t] mass obtained here has been 
compared with the previous measurements in fig. 
\6\ This result confirms the measured value of the 
NA48 at 0.6cr, with an error reduced by a factor 
2 and is lla away from the GEM result, it is 
also compatible, at 1.4cr level, with less accurate 
CLEO result [10]. 
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Figure 6. rj mass measurements, see text for the 
references. The continuous Une has been com- 
puted using the PDG procedure J^, pag. 14. 
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